Light emission patterns from stadium-shaped semiconductor microcavity lasers 
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We study light emission patterns from stadium-shaped semiconductor (GaAs) microcavity lasers 
theoretically and experimentally. Performing systematic wave calculations for passive cavity modes, 
we demonstrate that the averaging by low-loss modes, such as those realized in multi-mode lasing, 
generates an emission pattern in good agreement with the ray model's prediction. In addition, we 
show that the dependence of experimental far-field emission patterns on the aspect ratio of the 
stadium cavity is well reproduced by the ray model. 
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To theoretically investigate emission patterns from 
two-dimensional microcavity lasers, two approaches have 
been successfully employed, one based on the wave de- 
scription and the other based on the ray description. 
The ray description is regarded as an approximation of 
the wave description and can be justified in the short- 
wavelength limit by the Eikonal theory only in limited 
cases, i.e., when ray dynamics becomes integrable. It 
has been a fundamental problem in the field of quantum 
chaos to study how the nonintegrable or chaotic ray dy- 
namics manifests itself in the wave description [ij. In 
spite of the lack of full justification, the ray description 
has been practically used for cavities that exhibit chaotic 
ray dynamics, and shown useful for giving a simple ex- 
planation for the appearance of emission directionality 

Recent studies on microcavities with low refractive in- 
dices (e.g. n = 1.5 for polymers and n = 1.36 for dye- 
doped liquid jets) have revealed a remarkable ray-wave 
correspondence in low-loss cavity modes; the far-field 
emission pattern of an individual low-loss cavity mode 
exhibits close agreement with the result generated by a 
ray model [1, [a 0, B l^- So far, it is not clear whether 
this ray-wave correspondence is due to the low refrac- 
tive indices (i.e., high openness) or a more robust prop- 
erty holding also for higher refractive index cases. The 
present work examines the ray-wave correspondence for 
semiconductor (GaAs) cavities with the stadium shape 
[lo| . whose refractive index is much higher (n = 3.3) 
compared to polymers and dye-doped liquid jets. While 
there exist several experimental and theoretical works for 
GaAs microcavities [ll|, |l2|, |l3|, \l4 1 , systematic analysis 
has not yet been performed for low-loss cavity modes in 
the short-wavelength regime. 

In this paper, we show that the far-field emission pat- 
tern of an individual low-loss mode does not always ex- 
hibit good agreement with the ray model's prediction, 
being different from the cases for the low refractive index 
cases. Nevertheless, a new correspondence is found in the 
phase space distributions describing near-field emission 
patterns. We demonstrate that close correspondence be- 



tween the ray and the wave description in the far field is 
retrieved by performing the averaging by low-loss cavity 
modes. In addition, we show that an experimental far- 
field emission pattern, which can be approximated by the 
averaged result of low-loss modes, systematically agrees 
with the result of the ray model. 

We define the shape of a stadium cavity in the inset 
of Fig. [T] (b). We introduce an aspect ratio parameter 
e — L /R, where L is the half length of the linear part and 
R the radius of the circular part. We assume that the 
refractive index inside the cavity is n = 3.3 and n — 1.0 
outside the cavity. For e = 1.0, far-field emission patterns 
have been obtained experimentally and their trends have 
been explained by the ray model 



Whereas the far- 
field emission pattern is less directional for e = 1.0, by 
decreasing the value of e to less than around 0.3, one finds 
drastic changes in the far-field emission patterns. This 
drastic change can be associated with the phase space 
flow governed by the unstable manifolds of three-bounce 
periodic orbits, for which all the bounces occur at the 
circular parts. The details of this drastic change will be 
reported elsewhere. 

In Fig. [H we plot far-field emission patterns for e = 1.0 
and e = 0.3 generated by a ray model. In the ray 
model, we incorporate Fresnel's law for transverse mag- 
netic (TM) polarization to describe the light leakage at 
the cavity boundary [l5 , l^, 17 1 . Another method to 
theoretically obtain far-field emission patterns is based 
on wave calculation. Passive cavity modes are calculated 
by using, for instance, the boundary element method [lit . 
These are the eigensolutions of the Helmholtz equation 
[V^ -t- n^fc^]?/; = 0, where fc is a wave number outside the 
cavity and becomes complex by imposing the outgoing 
wave condition at infinity. For computational simplicity, 
we assume that the electric field is TM; we regard ip as 
the z component of the electric field and assume that ip 
and its normal derivative dtp are continuous at the cavity 
boundary. Because of computational limitation, we set 
the size parameter as kR « 100 (i.e., nkR « 330), which 
is about half of a real cavity size used in experiments 
discussed below. Since we have previously confirmed the 
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convergence of wave calculations even when nkR is less 
than 70 @], we think nkR « 330 is large enough to dis- 
cuss semiclassical properties and thus allows one to com- 
pare wave calculations with experiments. Also, we note 
nkR « 330 is much larger compared with recent previ- 
ous works 0, S El. For e = 1.0 and e 0.3, we ob- 
tained in each case around 100 passive cavity modes for 
99.95 < ReA:i? < 100.05. We are interested in the emis- 
sion patterns of low-loss modes, since it is these modes 
that contribute to lasing. We plot in Fig. [U far-field 
emission patterns for low-loss modes. For n = 1.5, it was 
relatively easy to confirm the correspondence between the 
results of wave calculations and those of the ray model, 
since the ray model predicts the appearance of four dis- 
tinct narrow peaks in the far-field emission pattern 0] ■ 
However, for the less directional cases shown in Figs. [T] 
(a) and[T](b), the correspondence between wave calcula- 
tions and ray simulations becomes less clear. In particu- 
lar, in Fig. [1] (a), we see a large discrepancy at 9 = 90°, 
where 6 is the far-field angle defined in the inset of Fig. 
[1] (b) . We note that the peaking at 6* = 90° is not a com- 
mon feature of low-loss modes. By averaging the far-field 
emission patterns over low-loss modes, we see that the av- 
eraged patterns correspond closely with the ray model's 
results. 

The correspondence between the ray model and the 
wave description can be further studied by looking at 
the near fields. Here we introduce functions describing 
near-field emission patterns, defined in the phase space 
spanned by the Birkhoff coordinates {s,p), where s is 
the curvilinear coordinate along the cavity boundary and 
p the momentum component tangential to the cavity 
boundary. 

In the ray model for stadium cavities, the light in- 
tensity inside the cavity £{t) can be asymptotically de- 
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scribed as £{t) oc exp[— 2 
7fl can be expressed as [l 

IR. = / ds 



l/r. 



where the decay rate 



l/n 



dp Pis, p), 



(1) 



where S is the total boundary length and we assume 
momentum is normalized to unity. The function P{s,p) 
describes how much light is transmitted outside the cav- 
ity at a boundary point s in the direction determined by 
p. We plot P{s,p) for a stadium cavity with e = 1.0 in 
Fig. [11(a). 

For a passive cavity mode, on the other hand, the cor- 
responding decay rate can be expressed as 



iw 



ds 



l/ri 



l/n 



dpH{s,p). 



(2) 



Here, H{s,p) is a Husimi-like phase space distribution 
calculated from the wave function at the cavity boundary 
^{s) and its normal derivative d^p{s), defined by 
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FIG. 1: (Color) Far-field emission patterns numerically cal- 
culated for n = 3.3 stadium cavities with (a) e — 1.0 and 
(b) e — 0.3. 6 is the far-field angle defined in the inset of 
Fig. [1] (b) . The blue curves are the patterns obtained from 
the ray model. The green curves are the patterns for cavity 
modes with low loss. The mode in (a) has the seventh low- 
est loss and the mode in (b) has the third lowest loss in the 
searched kR range. The red curves are the average of the 
far-field patterns of the 30 lowest-loss modes. We note that 
all the far-field emission patterns are normalized so that the 
integration becomes unity. 



where hf{s,p) = J ds'G* {s'; s,p)f{s') and G{s';s,p) is a 
coherent state for a onc-dimensional periodic system 



G{s';s,p) = 
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H{s,p) ^ Im[/i^,(s,p)/i9^(s,p)], 



(3) 



2^"^' — hip(s'-s-mS)| 

(4) 

with a = ^ S'/[2nRe(fci?)]. Below, we will investigate the 
correspondence between H{s,p) and P{s,p). 

In Fig. [11(b), we plot H{s,p) for a low-loss cavity mode 
whose far-field emission pattern is shown in Fig. [l] (a). 
We superimpose a set of points [red curve in Fig. [2 (b)] 
giving the far- field emission at = 90° . One can see that 
the strong near-field emissions indicated by arrows result 
in the strong far-field emission at 9 = 90°, which can be 
observed in Fig. [l] (a). 

In Fig. m (b), we can also see that H{s,p) is mainly 
supported on the high-intensity regions of P{s,p). This 
property has already been found for stadium cavities with 



3 




s / S 



FIG. 2: (Color) Phase space distributions for a stadium cavity 
with n — 3.3 and e = 1.0. (a) P{s,p) generated by the ray 
model, (b) H{s,p) for a low-loss cavity mode whose far-field 
emission pattern is plotted in Fig. [1] (a). The points on 
the red curve contribute to the far-field emission ai 9 = 90°. 
High-intensity regions resulting in strong far-field emission at 
9 = 90° are indicated by arrows, (c) The average of H{s,p) 
of the 30 lowest-loss modes. 



1.5 as a common feature of low-loss modes [17 1 



Below, we systematically check this property for n = 3.3 
by introducing a quantity measuring the overlap of the 
high-intensity regions of P(s,p) and those of 7J(s,p). We 
define the high-intensity regions of a function f{s,p) as 



where / is the average 
J^'dsj'_{';^dpfis,p)/i2S/n). 
and can be defined as 



of f{s,p), i.e., / = 
The overlap between Qp 



fi{np n fin) 
^i{VLp)ii{Q.h) ' 



(6) 



where ^{Vl) is the Lebesgue measure of the set 51. The 
overlap / takes the maximum value 1=1 when D,p = Qh 
and takes the minimum value 7 = when Qp n ^Ih = 0- 
For instance, the I value for H{s,p) shown in Fig. [5] (b) 
is 0.66. For e = 1.0 and e = 0.3, we plot the / values as 
a function of the loss rate |Im(fci?)| as shown in Fig. [31 
For the both e values, one can confirm a clear trend that 
the smaller the loss rate, the larger the overlap. 

As can be seen in Fig. [5| (b), H{s,p) is localized 
strongly on some portion of flp. This localization causes 
the discrepancy in far-field emission patterns between 
a cavity mode and the ray model. Since the localiza- 
tion pattern varies depending on the mode, by averaging 
H{s,p) over low-loss modes, one obtains a smeared dis- 
tribution H(s,p), which agrees more closely with P(s,p) 
than the individual H(s,p). This improvement of the 
ray-wave correspondence explains why the averaged far- 
field emission patterns shown in Fig. [T] correspond better 
with the results of the ray model. In Fig. [2] (c), we show 
the average of H{s,p) of the 30 lowest-loss modes. The 
overlap / between the high-intensity regions of H{s,p) 
and Op becomes 0.82 for e = 1.0 and 0.73 for e = 0.3, 
which are indicated by dashed lines in Fig. [3] 

Let us explain why it has been found for n = 1.5 that 
each of the low-loss modes always has a far-field emis- 
sion pattern closely corresponding to the results of the 
ray model d, 0I- In contrast to the case of n = 3.3, 
for n = 1.5 the high-intensity regions of P{s,p) turn out 
to consist of narrow stripes [l7|. Moreover, each of the 
stripes is parallel to the curves of the constant far-field 
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FIG. 3: (Color online) The overlap / of fip and ^Ih as a 
function of the loss rate |Im(fci?)| for e = 1.0 (left) and e = 
0.3 (right). The overlap of the high-intensity regions of the 
averaged H and f2p is indicated by a dashed line. 
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FIG. 4: (Color) Far-field emission patterns for stadium- 
shaped semiconductor microcavity lasers. Red and green 
curves are experimental data for two different samples fab- 
ricated in the same manner. Blue curves are data generated 
by the ray model. In all data, intensities are normalized so 
that the integration becomes unity. 



angle. As a result, irrespective of how H{s,p) is localized 
on rip, the localized regions are always located near the 
curves of the constant far-field angle, and thus generates 
a highly directional far-field emission pattern consistent 
with the result of the ray model. To summarize, several 
conditions need to be satisfied for the correspondence be- 



tween the far-field pattern of an individual low-loss mode 
and the ray description. Therefore, the ray-wave corre- 
spondence for an individual low-loss mode on the level 
of far-field emission patterns can be observed only for a 
specific choice of the parameters such as cavity shape and 
the refractive index. On the contrary, we expect that the 
agreement of the supports of P{s,p) and H{s,p) can be 
observed more robustly for an individual low-loss mode 
in a sufficiently semiclassical regime. 

So far, we have shown numerically that the ray model 
for n = 3.3 can be validated when one considers the av- 
erage of many low-loss cavity modes in the semiclassi- 
cal regime. We checked that this correspondence can 
be also observed between transverse electric (TE) wave 
calculations and the ray model with Fresnel' law for 
TE polarization, the detail of which is reported in Ref. 
[20| . Lastly, we examine whether the ray model can 
explain experimental far-field data for stadium-shaped 
GaAs microcavity lasers when multiple modes are in- 
volved in lasing. We fabricated stadium-shaped micro- 
cavity lasers by using a metalorganic chemical vapor 
deposition grown gradient-index, separate-confinement- 
heterostructure, single-quantum- well GaAs/Ala;Gai_a;As 
structure and a reactive-ion-etching technique ■ The 
radius R of the stadium is fixed as 25 ^m, and the e value 
is varied from 0.2 to 1.0. Lasing is achieved at room tem- 
perature by using a pulsed current with 500-ns width at 
a 1 kHz repetition. We confirmed the sharp narrowing of 
optical spectra around 850 nm above the lasing thresh- 
old, and we checked that lasing occurs in multimodes. 
For instance, for the cavity with e — 1.0 we found 13 
peaks in the spectrum. In Fig. [4] we plot measured far- 
field emission patterns. For each e value, we plot data for 
two different samples (red and green curves), which are 
fabricated in the same way and pumped with the same 
injection current strength. We also superimpose the far- 
field emission patterns obtained from the ray model in the 
blue curve. In the ray model we employed Frenel's law for 
TE polarization, taking into account that experimentally 
observed emission is TE polarized. The resulting TE far- 
field patterns of the ray model are slightly different from 
those of TM polarization. In Fig. [H one can confirm 
that changes in the experimental far-field emission pat- 
terns due to the change in e value are well explained by 
the ray model. 

In general, the relation between cavity modes and 
a multi-mode lasing state formed by nonlinear mode- 
interaction through an active medium can become very 
complicated However, the present work revealed 

that, at least concerning far-field emission patterns, the 
average of low-loss modes approximates a multi-mode las- 
ing state for a sufficiently large nkR value. We believe 
that the enhancement of the ray-wave correspondence 
due to the averaging by low-loss modes plays a signif- 
icant role in yielding the good correspondence between 
the experimental far-field emission patterns and those of 
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the ray model. 

In summary, we studied light emission patterns of low- 
loss modes for stadium-shaped semiconductor microcav- 
ities. The correspondence between low-loss modes and 
the ray model was revealed by investigating the phase 
space distributions describing near-field emission pat- 
terns. Close correspondence was found between exper- 
imental far-field data and the results of the ray model, 
which we attributed to the enhancement of the ray- wave 
correspondence due to the averaging by low-loss modes. 

The work at ATR was supported in part by the 
National Institute of Information and Communications 
Technology of Japan. 
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